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One of the solutions to the cosmological Polonyi problem is to introduce a large 
coupling between the Polonyi field and the inflaton so that the Polonyi field adiabat- 
ically tracks the temporal minimum of the potential. We study general conditions 
^ i for the adiabatic suppression mechanism to work, and find that a non-negligible 

^ \ amount of the Polonyi field is induced in the form of coherent oscillations at the 

end of inflation. In the case of low reheating temperature, this contribution is so 
small that it does not cause cosmological problems. On the other hand, this con- 
tribution may be significant for a relatively high reheating temperature and we still 
need some amount of tuning in order to avoid the Polonyi problem. We also point 
out that Polonyi particles produced from thermal plasma pose a severe constraint 
on the reheating temperature. Furthermore, we extend the original framework to 
include enhanced couplings of the Polonyi field with the visible particles as well as 
with itself, and derive upper bounds on the reheating temperature after inflation. 
We also investigate the adiabatic solution to the cosmological moduli problem in 
gauge and anomaly mediation. 



1 Introduction 



Supersymmetric (SUSY) extension of the standard model (SSM) is a plausible candidate 
for physics beyond the standard model since it solves the gauge hierarchy problem and 
leads to the successful gauge coupling unification at the grand unification scale. It also 
contains candidates for dark matter (DM). 

The gravity-mediation models for SUSY breaking are attractive because of its simplic- 
ity: the SSM gaugino masses arise from an F-term of an elementary gauge singlet field Z 
through Planck-scale suppressed interactions. However, the Polonyi field Z is known to 
cause a serious cosmological problem. Since it is neutral under any symmetry, there is no 
special point in its field space. Therefore, the minimum of the effective potential for the 
Polonyi field during inflation is generically deviated from the one in the low energy. After 
inflation, the Polonyi starts to oscillate about the minimum with an amplitude of order 
the Planck scale Mp, and soon dominates the energy density of the Universe. Since its 
interactions are suppressed by the Planck scale, the lifetime of the Polonyi is very long, 
leading to an onset of a radiation-dominated Universe with a low temperature, typically 
below MeV. Such a low temperature would dramatically alter the big bang nucleosyn- 
thesis (BBN) predictions of light element abundances in contradiction with observations. 
This is the notorious cosmological Polonyi problem [T]. 

Even if one gives up the gravity mediation, there may be scalars having similar prop- 
erties to the Polonyi field. Indeed, there generally appear such scalars, called moduli, 
associated with the compactification of extra dimensions in the string theory. The moduli 
have a cosmology similar to the Polonyi, and cause a serious cosmological problem [2]. 

Some of the solutions to the cosmological Polonyi/moduli problem require significant 
modification of the conventional cosmological scenarios. One of them is to make the 
Polonyi/moduli heavy enough to decay well before the BBN begins. In the anomaly- 
mediated SUSY breaking models the gravitino as well as the moduli are expected 
to have masses much heavier than the weak scale. However, the solution turned out to 
be more complicated, because of unsuppressed production of SUSY particles including 
gravitinos from the modulus decay P]. The modulus decay also significantly dilutes the 
pre-existing baryon asymmetry. This makes the most baryogenesis mechanisms unsuc- 
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cessful, although it is still possible to create the right amount of baryon asymmetry by the 
Affleck-Dine mechanism [5l[6]. Another solution is to introduce late-time entropy produc- 
tion for diluting the Polonyi/moduli abundance. For instance, the thermal inflation [7] 
can provide enough dilution. Again, it also dilutes any pre-existing baryon asymmetry to 
a negligible amount. Considering that thermal inflation should dilute the modulus abun- 
dance by up to a factor of ~ 10^° in order to solve the moduli problem, it is hopeless to 
try to create a sufficient amount of baryon asymmetry which survives the huge dilution. 
Thus we need to create the baryon asymmetry after thermal inflation by an elaborate 
mechanism [8]. 

An interesting solution to the Polonyi/moduli problem was proposed long ago by 
Linde [9]. It was pointed out that, if the modulus fleld has a mass squared of ~ (?H'^ before 
it starts to oscillate, with c = (9(10) and H being the Hubble parameter, the modulus 
follows the time-dependent potential minimum adiabatically and the resultant amplitude 
of coherent oscillations is signiflcantly suppressed. This simple solution is attractive since 
it works for a wide range of the Polonyi/modulus mass, and since there is no need for 
an additional late-time entropy production, which may make the standard leptogenesis 
scenario [10] viable. Recently, two of the present authors (FT and TTY) noticed that there 
is an upper bound on the reheating temperature for the adiabatic solution to work [11] and 
also showed that such a large Hubble mass may be a consequence of the strong dynamics 
at the Planck scale [12] or the fundamental cut-off scale one order of magnitude lower 
than the Planck scale [TT] . 

In this paper we study the adiabatic suppression mechanism in a great detail in or- 
der to establish the solution in a complete form and explore parameter space where the 
Polonyi/moduli problem is solved. First, we carefully investigate the issue of adiabatic- 
ity and how the Polonyi/moduli oscillation is induced after inflation. We flnd that the 
coherent oscillations of the Polonyi/modulus fleld are generically induced at the end of 
inflation where the adiabaticity is violated and the produced amount depends on infla- 
tion models. The adiabatic solution may not work without tuning parameters or an 
additional entropy production for a relatively high reheating temperature, although the 
required amount of tuning or entropy production is greatly relaxed with respect to the 
original Polonyi/moduli problem. For a low reheating temperature, on the other hand. 
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the adiabatic solution can solve the Polonyi/moduli problem without any fine-tuning. 

Next, we point out that, even if coherent oscillations of the Polonyi/moduli are 
reduced to a negligible amount by the adiabatic solution (and some fine-tuning), the 
Polonyi/moduli are generically produced by particle scatterings in thermal plasma as 
long as they have (Planck-suppressed) interactions with the SSM particles. In particular, 
the thermal production of the Polonyi field is inevitable because it must be coupled to the 
visible sector to mediate the SUSY breaking. We find that the thermal production of the 
Polonyi/moduli fields leads to non-trivial cosmological constraints. In the gravity media- 
tion, the constraints are so stringent that there is no parameter space where the thermal 
leptogenesis scenario works, other than the heavy gravitino mass region, m^i2 ^ lOTeV. 

Then, based on the findings of Refs. [121 E]) we extend the original framework to 
include an enhanced coupling of the Polonyi field with itself as well as the SSM particles. 
We will see that a new interesting possibility emerges in this case, where the gravitino 
is the lightest SUSY particle (LSP) and relatively high-reheating temperature is allowed. 
The focus-point like mass spectrum is favored in this case. 

We will also discuss the implications of the adiabatic solution to the moduli problem 
in the gauge-mediated SUSY breaking models where the moduli and gravitino are light. 
For light moduli, < 1 MeV, the adiabatic suppression mechanism can solve the mod- 
uli problem without fine-tuning on the modulus potential. However, we show that the 
reheating temperature is severely bounded from above due to the gravitino thermal pro- 
duction and the condition for the adiabaticity, and the thermal/non-thermal leptogenesis 
does not work even in a scenario with an extremely light gravitino (^3/2 ^ 16 eV |13])F1 

Lastly we will discuss the adiabatic solution to the moduli problem in the anomaly- 
mediated SUSY breaking models where the moduli and gravitino are heavy enough to 
decay before BBN. In this case, the adiabatic suppression can solve the moduli problem 
for a relatively high reheating temperature once we allow a tuning of (9(0.01) on the 
modulus potential. The (non-)thermal leptogenesis may work for the gravitino mass of 
~ 100 TeV with the Wino LSP with mass of a few hundred GeV. 

^ The extremely light modulus with a mass smaller than 0(1) keV in the form of coherent oscillations 
contributes to the DM density of the Universe. Hence anthropic arguments may guarantee the smallness 
of the modulus abundance. 
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The rest of this paper is organized as follows. In Sec. [2] we investigate the adiabatic 
suppression mechanism, and estimate the abundance of the Polonyi/moduli field in the 
form of coherent oscillations. We also derive an upper bound on the reheating temperature 
for the mechanism to work [11]. In Sec. |3]we study the Polonyi/moduli oscillation induced 
at the end of inflation in detail and show that it significantly contributes to the final 
Polonyi/moduli abundance. We derive the required amount of fine-tuning to solve the 
moduli problem. In SecHlwe study the Polonyi problem in the gravity mediation in detail 
in order to clarify to what extent the Polonyi abundance is suppressed in the mechanism, 
taking account of thermal production of the Polonyi field. In Sec. [5] and [6] we will discuss 
the moduli problem in gauge and anomaly mediation. The last section is devoted to the 
discussion and conclusions. 

2 Adiabatic solution to the Polonyi/moduli problem 
2.1 Basic idea 

First we briefly review the basic idea to suppress the modulus abundance Let z 
denote collectively a modulus field, including the Polonyi field. We set the origin of z 
so that it coincides with the potential minimum at present. In the early Universe, the 
effective potential of z could receive various corrections from its interactions with the 
infiaton and/or SSM particles. In particular, if z has a quartic coupling with the infiaton 
in the Kahler potential, z receives the so-called Hubble-induced mass term. Here and in 
what follows, the infiaton also refers to a field which dominates the energy density of the 
Universe when the modulus starts oscillating. 

Let us express the effective potential in the early Universe as 

V = \my + \c'H\z - z.)\ (1) 

where is the mass of z in the low energy, z^, is the initial displacement during inflation, 
H is the Hubble parameter and c is a constant. The amplitude 2;* is expected to be of 
order the Planck scale. For simplicity we treat 2; as a real scalar field when we consider 
its dynamics, but this does not affect the following argument. We have assumed that the 
potential can be approximated by a quadratic potential about the origin, at least up to 
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z = z^. This is expected to be the case for the Polonyi field, but the potential may take 
a more general form. We will come back to this issue in Sec. [31 

Let us first consider the case of c = 0{1). This is the case if the modulus and the 
inflaton are coupled by the Planck-suppressed operators with coefficients of order unity. 
The modulus dynamics is as follows. Assuming that the Hubble parameter during inflation 
is much larger than m^, z is stabilized at 2; ~ during inflation. After inflation, it begins 
to oscillate around the minimum z = with an amplitude of when H ~ rriz- The 
modulus abundance in the form of coherent oscillations is given by 



where pz is the modulus energy density, s is the entropy density, r$ is the inflaton decay 



Tosc = {lO/n'^g^y^^VmzMp. The above modulus abundance is so large that it causes a 
serious cosmological problem. For example, we need a tuning of 2* < lO^^^Mp to satisfy 
the BBN bound for a typical reheating temperature Tr = 10^ GeV and rriz ~ 1 TeV. If c 
is much smaller than 0{1), the modulus abundance depends on the initial displacement, 
which is subject to quantum fluctuations during inflation. As long as it is of order z^,, the 
resultant modulus abundance is the same order as in the case of c = 0{1). 

The situation significantly changes if c ^ One might expect that the modulus 

begins to oscillate when H ~ c^^rriz with an amplitude of z*, when the potential minimum 
starts to move from to the origin. However, this is not true. At i7 ~ c~^m^, the Hubble 
parameter is much smaller than rriz, which means that the potential minimum moves more 
slowly than the typical time scale for the modulus dynamics. Thus there is enough time 
for z to follow the potential minimum, and as a result, the coherent oscillations of z are 
significantly suppressed. 

It is possible to estimate the suppression factor analytically, by solving the equation 
of motion for the modulus in the potential ([1]). The final amplitude of the coherent 
oscillations is suppressed by the following factor [9j 




for < TUz 





rate, Tr = ^JT^Mp denotes the reheating temperature after infiation, and 
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where p parametrizes the Hubble parameter as H = p/t when H = c~^mz, and it is given 
by p = 2/3 and p = 1/2 before and after the reheating, respectively. Thus, as we increase 
c for fixed and 2;*, the final modulus abundance is exponentially suppressed by 
relative to the estimate (E]). In Fig. [T]we have shown the numerically obtained evolution 
of the modulus abundance as a function of time as well as the analytic estimates based on 
([2]) and ([3]) for c = 1 (top) and 10 (bottom). Here we took z^, = Mp/c. We can see that 
in both cases the analytic estimates agree well with the numerical results. In this plot, 
the entropy density s is defined as s(t) = s(TR)[a{T^) / a(t)]^ and s(Tr) = (27r^(yf*s/45)T^. 
The definition of s, and hence Pz/s, coincides with the standard one at late times after 
reheating. 
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Figure 1: Time evolution of the modulus abundance for c = 1 and c = 10. For comparison. 



analytic estimations, labeled as (ana), are also plotted. We have taken 
and = Mp/c. 
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Thus, we can achieve sufficient suppression of the modulus abundance for c = 0(10) 
because of the exponential factor appearing in Eq. ([3]), which reflects the fact that the 
variation of the adiabatic invariant is exponentially suppressed. This is essentially what 
was pointed out in the pioneering paper [9j. 

Lastly, we note here that the above estimate ([3]) is based on the analysis of Ref. [H], in 
which the initial condition is set to be z = and i = at t = when the Hubble param- 
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eter is infinitely large. However, we would like to emphasize here that the actual initial 
condition of the modulus field is set during inflation when the Hubble parameter is finite. 
In addition, the infiaton dynamics changes the modulus potential at the end of inflation 
with a time scale of the infiaton mass, which may violate the adiabatic condition. There- 
fore, the above analytic estimate is only approximate, and there is in general additional 
contribution to the modulus abundance produced at the end of inflation. That is to say, it 
is important to follow the modulus dynamics during the first oscillation At = 0{{cH)~^) 
after inflation, and non-negligible amount of the modulus abundance can be generated 
during the short period. Once the modulus starts to oscillate, the particle production is 
indeed exponentially suppressed as explained above. This additional particle production 
at the end of inflation has been overlooked so far in the context of the adiabatic sup- 
pression mechanism, although similar kind of particle production was considered in the 
detailed analysis of thermal inflation [H]. We will come back to this issue in Sec. El and 
derive the condition for the adiabatic suppression to work successfully. For the moment, 
we will neglect the additional contribution. 

2.2 The origin of the enhanced coupUngs 

The essential ingredient for the mechanism described above is the large Hubble mass term 
with c ^ 1. Suppose that the infiaton <I> has a quartic coupling with the modulus in the 
Kahler potential as [9] 

K ^ -c?^. (4) 

where Z and $ are chiral superfields corresponding to the modulus and the infiaton, 
respectively. The c is related to ci as 

c2 = 3(c2 + l). (5) 

If Ci is larger than order unity for some reason, a large Hubble-induced mass term for 
the modulus is generated. The origin of Ci ^ 1 may be due to a strong dynamics at the 
Planck scale p[2] or a cut-off scale one order of magnitude below the Planck scale 

We would like to emphasize here that the Hubble-induced mass term generated by 
(jl]) disappears after the reheating, since the energy density of the Universe becomes 
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dominated by radiation afterward. Therefore, the reheating should occur sufficiently 
late for the adiabatic solution to work. This places an interesting upper bound on the 
reheating temperature [11], which we shall derive explicitly in the next subsection. 

We note here some subtleties regarding the thermal effects. If the modulus is coupled to 
the SSM particles in plasma, the modulus potential receives thermal corrections after the 
reheating. In particular, if the thermal mass is much larger than the Hubble parameter, 
the modulus adiabatically follows the time-dependent minimum of the potential in a 
similar way discussed above. The thermal mass is roughly estimated as 



ml ~ = 0(0.01) X d^H\ (6) 

where we have used H"^ = {n'^g^/90)T'^/Mp with ~ 230 in the second equality. Here c' is 
the coupling constant between the modulus and SSM particles in thermal bath normalized 
by the Planck scale; it corresponds to e.g. C3 or C5 in Eqs. f l53|) and fl5^ . Thus, if c' is larger 
than (9(100), the adiabatic solution may work even if the reheating is completed before 
H = rriz/c. However, as we shall see later, such a large coupling between the modulus and 
the SSM particles leads to problematic thermal production of the modulus. Therefore we 
assume c' < (9(10) in the following so that the thermal mass is at most comparable to the 
Hubble mass with c = 0{1) and it does not affect the modulus dynamics significantly^ 

2.3 Upper bound on the reheating temperature 

Assuming that the large Hubble-induced mass term arises solely from an enhanced cou- 
pling between the modulus and the infiaton, we can derive an upper bound on the re- 
heating temperature for the solution to work. It is important to note that the large 
Hubble-induced mass term disappears as e"^"**, where denotes the decay rate of the 
infiaton. The temporal minimum moves on a time scale of F^^ at the reheating, which 
should be smaller than for the adiabatic solution to work. Therefore, we have [TT] 

< m,. (7) 

It depends on the required suppression factor as well as on the value of c how much the 
decay rate should be suppressed compared to the modulus mass. From Fig. [21 one can see 



^ There is also a linear term such a.s V ^ T*z/M [T5], whose effects on the modulus abundance were 
discussed in Ref. [IB] . 
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that the adiabatic solution does not work unless < O.lm^, and that the effect of the 
reheating becomes practically negligible if < O.Olm^. We therefore adopt < 0.05 
as a reference value in the following analysis. 

We can express the bound in terms of the reheating temperature Tr as [H] 



Tr < 3xl09GeVf^r^f^^)\ (8) 
^ ~ V230/ VlTeV/ ^ ^ 

where the reheating temperature Tr is related to the decay rate as 

Tn^(^^y"Vm^- (9) 

The upper bound on the reheating temperature was overlooked so far, but it has a very 
important implication especially for the baryogenesis scenario such as thermal leptogen- 
esis. 

2.4 Suppressing the modulus abundance 

We numerically evaluate the suppression factor for the modulus abundance, A, which is 
defined by the ratio of the actual modulus abundance to the analytic estimate for the case 
of c = 1, 

\Pz/ ''/cxp 

where the denominator is given by Eq. ([2]). 

Contours of the suppression factor A in the plane of (c, F^/m^) are shown in Fig. |2j 
The initial amplitude is taken to be z^, = Mp/cjj It is seen that the modulus abundance 
is highly suppressed for c > 30 and T^/niz < 0.05. Practically, A < 10^^° is sufficient 
for satisfying the BBN bound on the modulus abundance, that is, pz/s < 10~^^GeV 
for rriz ~ 1 TeV [T7]. With this amount of suppression, the moduli do not dominate 
the Universe and therefore do not produce huge entropy with their decays. Thus the 
baryon asymmetry is not diluted. On the other hand, if we merely demanded that the 
moduli do not cause entropy production, the required suppression factor would be mild 



The large c effectively corresponds to a small cutoff scale of order of Mp /c in the non-renormalizable 
Kahler potential. Thus the amplitude should be smaller than or comparable to this cutoff scale in order 
for the effective description to be valid. 
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(say, A < 10"^° for = ITeV and Tr ~ lO^GeV). Hereafter we regard A ~ lO^^o 
as a typically required suppression factor for solving the cosmological moduli problem, 
although the precise constraint depends on the modulus mass and couplings, and also the 
reheating temperature. 
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Figure 2: Contours of the suppression factor A in the plane of (c, F^/m^). 



3 When does (not) the adiabatic solution work ? 

We have seen that the adiabatic solution works for c = (9(10) in the potential (II]). How- 
ever, there is in general an additional contribution to the Polonyi/modulus abundance, 
and it is necessary to clarify under which conditions the adiabatic solution works. 

Before going to concrete models, we give a general discussion on the adiabatic sup- 
pression mechanism. Suppose that the scalar potential for a modulus z depends on the 
Hubble parameter. During inflation the z is stabilized at a point deviated from the 
low-energy minimum, and the position of the minimum gradually changes as the Hubble 
parameter decreases after inflation. Then question is in which situation the coherent os- 
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cillations are induced. The key concept is the adiabaticity of the modulus dynamics in 
the time- dependent scalar potential. Let us assume that a modulus with a mass larger 
than the Hubble parameter is oscillating about the potential minimum z^[^{t). If the rate 
of the change of the potential minimum is much smaller than the oscillation frequency, 
the modulus number density in a comoving volume is conserved. The condition is written 
as 



where Sz denotes an amplitude of oscillations about the potential minimum. If the poten- 
tial can be approximated with a quadratic potential about Zmin, the far right hand side 
becomes equal to the mass. If this condition is violated, particle production occurs as 
2;min movcs. Thus the condition ( ITTl) must be satisfied in order for the modulus amplitude 
to be suppressed. 

A couple of comments are in order. First, the adiabatic suppression mechanism relies 
on the conservation of the adiabatic invariant (~ the modulus number density in the 
comoving volume), which is defined for a periodic motion. In other words, it does not say 
anything about the initial abundance before the Polonyi/modulus field starts oscillating, 
which is determined by the dynamics during the first period of oscillation. Therefore, 
careful case-by-case analysis is needed in order to estimate the initial Polonyi/modulus 
abundance. 

Secondly, if imin suddenly changes, the above condition ( ITT]) is not sufficient. (We here 
assume that 2;min does not jump at the end of inflation.) In this case, one should also 
consider a condition similar to f lTTj) with f{H) replaced with the typical time scale of the 
dynamics which determines the evolution of the potential minimum. For instance, as we 
shall see below, in a broad class of inflation models, the potential minimum starts to move 
at the end of inflation in a time scale of the infiaton mass. In this case one should regard 
the infiaton mass as the typical time scale. Then the adiabaticity condition is usually 
violated at the end of inflation, and coherent oscillations of the Polonyi/modulus field are 
induced. We shall study this effect in Sec. 13.11 

Another issue to be addressed is whether z^ain remains same after inflation. As we shall 
see in Sec. 13. 2[ in realistic single-fleld inflation models in supregravity, because of the non- 




(11) 



12 



trivial couplings between the inflaton and the modulus in the supregravity potential, the 
position of the minimum generically changes in a non-adiabatic way, and some amount 
of modulus oscillation is induced at the end of inflation. We will also see in Sec. 13.31 that 
in a class of multi-field inflation models, one (or more) of the fields acquires a non-zero 
F-term at the end of inflation, which in general changes the position of the potential 
minimum. For instance, in the hybrid inflation, the waterfall field has a vanishing F-term 
during inflation, but it starts to oscillate after inflation, acquiring a sizable F-term. The 
typical time scale of this phenomenon is given by the mass of the waterfall field. If this 
is shorter than the modulus oscillation period at that time, the adiabaticity condition is 
violated, unless the modulus is coupled to the waterfall field exactly in the same way as 
to the inflaton field. Thus, the adiabatic suppression mechanism becomes inefficient and 
the modulus abundance is not exponentially suppressed. 
In the following we study these issues in detail. 

3.1 Polony i /moduli production at the end of inflation 

Let us first consider the simple model ([1]), with the initial condition given at a finite 
Hubble parameter. Hi = p/ti. The essential difference from ([3]) is that, after inflation, the 
potential minimum shifts from by a finite amount during the first period of oscillation 
~ {cHi)~^. Note that the shift would be infinitely small if the initial condition was 
given at t = when the Hubble parameter is infinitely large (see Eq. ( IT4l) ). Once the 
modulus starts to oscillate about the time- dependent minimum, the particle production is 
suppressed exponentially afterwards. So let us study the dynamics during the first period 
of oscillation after the inflation. 

Suppose that the minimum changes by 6z within 6t = {cHi)~'^ after inflation. If the 
potential minimum changes proportionally to time in the first period of oscillation, the 
initial amplitude is given by 



where 2;min(^) denotes the time-dependent potential minimum, and the dot means the 
derivative with respect to time. The initial abundance of the modulus field created at the 
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end of inflation is therefore given by 

^min 

{u)y f 6z{ti 



8 \ Mp J V^min(t 



(13) 



In the model ([T]), we have 



where we used rriz ^ cifj in the second equalities. The inflaton-matter dominated phase 
is assumed in Eqs. f lT5]) and ffT^ . (See discussion in Sec. 12.31 ) The initial abundance for 
the model ([1]) is given by 

Note that the condition (fTTj) is satisfied if c ^ 1, and therefore the particle production 
afterward is exponentially suppressed. The final abundance is given by the sum of the 
estimate given in Sec. 12.11 and f|T6|) . Compared to {pz/s)cxp given by Eq. ([2]), there is a 
suppression factor, {m-z/cHiY in (fT6|) . Thus, the modulus abundance is power suppressed, 
if this additional contribution is dominant. 

We have checked that this analytic estimate agrees well with the numerical results. 
In Fig. [3l the dependence of the modulus energy abundance on the initial condition Hi 
for c = 20 (top) and c = 30 (bottom) are shown. Numerical solutions are compared with 
analytic estimate of the modulus abundance based on ([3]) labeled as "analyticl", and that 
of the initial modulus abundance based on (fT6|) labeled as "analytic2". We have taken 

= 10~^'^m2. It is seen that analytic estimates well reproduce the numerical results 
and that the initial modulus abundance flT6|) becomes non- negligible as the ratio Hi/m 
decreases. This validates our consideration above. 

In order for the above modulus abundance to be consistent with observations, the 
Hubble parameter during (precisely speaking, at the end of) inflation, ifjnf ^ Hi, should 
satisfy 

«..,>500™,(|;)-Y^) (17) 
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Figure 3: The dependence of the modulus abundance on the initial condition Hi for c = 20 
(top) and c = 30 (bottom). Numerical solutions are compared with analytic estimate of 
the modulus abundance ([3]) labeled as "analytici" , and that of initial modulus abundance 
f|T6l) labeled as "analytic2". We have taken = lO'^-^m^. 
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where A denotes the required suppression factor for the modulus abundance to be con- 
sistent with observation, normahzed by {pz/s)exp in Eq- Although this inequality 
is satisfied for many inflation models, it certainly places a non-trivial lower bound on 
the inflation scale. The conditions (ITTll and ([8]), together with c ^ 1, are the necessary 
conditions for the adiabatic solution to work in the model ([!]). 

In deriving (fT6!) . we have implicitly assumed that the modulus does not "see" the 
infiaton oscillations. This is a valid assumption for low-scale inflation models such as new 
inflation models, in which ^ Hi^f is satisfied where is the infiaton mass around the 
minimum. In this class of models, the modulus minimum starts to move soon after the 
inflation ends, and the displacement during the first modulus oscillation determines the 
initial abundance. The particle production at the end of inflation can be understood be- 
cause there is another time scale, i.e., the infiaton mass, which is lighter than the modulus 
mass during inflation but becomes much heavier after the end of inflation. As a result, the 
minimum to start moving in a time scale of m<^. In other words, the decoupling process 
of the infiaton violates the adiabaticity condition, leading to the particle production. 

The situation is slightly different for the chaotic inflation model where ~ Hi^{ and 
hence ^ ci^inf for c ^ 1. In this case the modulus remains heavier than the infiaton 
for a while after inflation, and the adiabaticity is violated at ~ /fad = ^(p/c. The 
modulus abundance produced at i^ad is estimated in a similar manner by replacing 
with ifad, and the result is 



We have confirmed that this agrees with numerical calculation. See Fig. |H Here the sup- 
pression factor is given by dlrriz/m^)^, with respect to (E]). Recalling that the suppression 
factor of ~ 10~^° is needed to solve the moduli problem, and that the infiaton mass in the 
single- field chaotic inflation model is given by ~ 2 x 10^^ GeV, the suppression factor 
(m^/m^)^ is so small that the moduli produced in this way is cosmologically harmless. 

3.2 Polony i/moduli production in single-field inflation models 

So far we have assumed that the position of the minimum does not jump at the end of 
inflation. In supergravity based inflation models, however, this is not necessarily the case. 




(18) 
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10^ 

Figure 4: The modulus abundance in the single-field chaotic infiation model. We have 
taken c = 20, = O.lMp, = 10~^Mp, = lO'^'^m^. Together shown is the 
analytic estimate based on ( IT8|l . 

We show below that the minimum actually changes at the end of infiation in general 
single- field infiation models in supregravity, including the model of Ref. [18]. 
For illustration, let us consider the Kahler potential as 

K=\Z\^ + \4>\' + c'\^^^^, (19) 

where c ^ 1 and (p denotes the infiaton superfield. We generically obtain the following 
form of the scalar potential from this Kahler potential, 

V = 3Hl\z\'' + c^H^\z-z,\^ (20) 

where 3H^M^ = with being the infiaton potential energy, and 3H^Mp = |0p + V^. 
The point is that the first term in (120|) depends only on the infiaton potential energy and it 
rapidly oscillates after infiation. The time averaged value of is given by {H'^) = 
if the infiaton behaves as matter during the oscillation around its minimum. Therefore, 
the position of the minimum during infiation is given by z^[^ = (?z^l[(? + 3), while the 
minimum after infiation is given by z^j^ = (?z^j{(? + 3/2). Since the time scale of the 
oscillation of the if| term, which is determined by the infiaton mass, is much greater 
than the modulus mass, the modulus does not feel the rapid oscillation of if? itself but 
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its time averaged value. Thus the modulus oscillation is induced just after inflation with 
an amplitude of 

The resulting modulus abundance is estimated as 

(») n / \ 2 / 

Pz 9 ( \ ( m 



Compared with the standard estimation (|2]), the abundance is suppressed by the factor 
~ mz/{c^Hi^{). This suppression is not sufficient, although the amount of tuning required 
for solving the moduli problem is significantly relaxed. For example, for c = 30, Tr = 
lO^GeV, m, = 1 TeV and Hi^f = lO^^ QeV, the tuning of \z^\ < IQ-^Mp is needed for 
solving the moduli problem. The estimate (!22|) is similar to that in the multi-field infiation 
model discussed below. We will show the validity of this estimate numerically in the next 
subsection. 

3.3 Polony i/moduli production in multi-field infiation models 

Now let us discuss multi-field infiation models, where the position of the minimum could 
drastically change at the end of infiation. The time scale of the change is usually de- 
termined by the mass of fields acquiring the F-term after infiation, e.g., the mass of the 
waterfall fields in the hybrid infiation. Let us consider a class of infiation models with a 
superpotential 

W = Xf{(P). (23) 

Most SUSY infiation models including hybrid infiation [T3|, smooth- hybrid infiation [20 , 
two-field new infiation [2T] chaotic infiation [22] and its variants [23] fall into this categoryjj 
Here it is the X whose F-term gives the infiaton potential energy during infiation. In 
the two-field new and chaotic infiation models, (p is regarded as the slowly-rolling infiaton 
field, while X does not participate in the infiaton dynamics. In the hybrid infiation model, 
X plays the role of the slowly-rolling infiaton, while is the waterfall field. 
In general, X and can have different couplings to the modulus field. 



T^-~2 \X\^\Z-Zx\^ . \(P\^\Z - Z^\^ 

^"''^ Ml ^'''^ Ml ' ^^^^ 



This form of the superpotential has been recently studied in detail in Ref. 
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in the Kahler potential. This yields the scalar potential for the modulus as[f| 

V=^[^x{\X\' + \Fx\') \Z -Z^f + cl + \Z - Z,f] , (25) 

where X and Z represent the lowest component of the corresponding superfields. During 
inflation it is Fx that contributes to the Hubble parameter, while and the kinetic 
energy of X and are subdominant or negligibly small. In the two-field new and chaotic 
inflation models, X sits at the origin and does not participate in the dynamics. After the 
inflation ends, both |0p and oscillate with the same amplitude and opposite phase, 

while |Xp and {F^l"^ remain negligible. On the other hand, in the hybrid inflation model, 
all of them oscillates with a time scale of the inflaton mass in a complicated manner. 

In order to take into such behavior in the multi-field inflation models, let us parametrize 
the modulus potential as 

V = \my + \^xHl{z - zxf + \clHl{z - z,)\ (26) 

where the Hubble parameter is divided into two parts as iP = Hf + We assume that 
Hf ~ and if| ~ during inflation and that Hf and if| oscillate with an opposite 
phase after inflation. Note that their sum is equal to the Hubble parameter and so it does 
not oscillate. The time-averaged values of Hf and i7| are given by 

(H^,) = {Hi) = ^. (27) 

During inflation, where the kinetic energy of the inflaton is negligible, the modulus sits 
at z = zx- After inflation ends, the temporal minimum oscillates around between zx and 
Z(jj, and it is given by 

CxHfzx + clHlz^ 



Its time derivative is 



-mm 2 tt2 , 2 tt2 



(28) 



[zx-z^), (29) 



(4ff^ + c^/j|)2 

where we have approximated m\ <^ c\Hf + c^H^. Hereafter we consider the case of 
cx = 0(10) and cx ^ since otherwise it is obvious that the coherent oscillations of 



^ There are other terms for the modulus potential, but is dominant if , ^ 1 and if the other 
terms are suppressed by the Planck scale with coefficients of order unity. 
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the modulus with an amphtude of ~ zx or is induced at the end of inflation, and the 
cosmological moduli problem is not solved if zx ^ z^ ^ Mp. 

Let us first consider the case of cxHi^f ^ m^, which is satisfied in the hybrid and 
new inflation models. In this case the modulus mass is much smaller than the typical 
frequency of the inflaton oscillations, and so it only feels the time-averaged potential. At 
the end of inflation, the potential minimum changes instantly (in a time scale of m'^^) 
from Zx to z^m, where the time-averaged potential minimum, z^i^ is given by 

Zmin ^ ic^xzx + clz^)/ic^x + cl). (30) 

Here we have used Eqs. ( |271) and ( l28ll . Thus the initial amplitude of the modulus oscilla- 
tion is 

5z{ti) r-u Zx- Z^in |-(^X - Z^)- (31) 

The modulus abundance produced at the end of inflation in this case is estimated as 



V - K-^) KM,) ■ '''' 

Although this may be still too large to avoid the cosmological Polonyi/moduli problem, 
the abundance is suppressed by a huge factor with respect to the original case without 
the adiabatic suppression mechanism. For example, for cx = 30, c^f, = 1, Tr = 10^ GeV, 
TJiz = 1 TeV and Hi^f = 10^^ GeV, the tuning of \zx — z^\ ~ lO^^Mp is sufficient for 
solving the moduli problem. 

In order to check this estimate, we have evaluated the modulus abundance in the 
SUSY hybrid inflation model with the superpotential 



W = KX{(j)(j)- li'^). (33) 
For simplicity we assume the following form in the potential (126|) . 



^1 - ^772 ' l34) 



3M|, 



rr2 /.^|Xp(|0p+|0p) + |0p + ,,, 

= 3M| • 



The modulus abundance in this model is shown in Fig. 13 We have taken k = 1, fi = 
O.lMp, = lO-^Mp, = 10^2.5^^ ^^^^ ^ Mp/cx and z^ = zx/2. Results for 
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(cx, c^) = (30, 2) and (40, 5) are shown together with the analytic estimate based on (!32ll . 
We can see that the analytic estimate on the final modulus abundances agree well with 
the numerical results. 




Figure 5: The modulus abundance in the hybrid inflation model. We have taken k = 1, 
fi = O.lMp, = IQ-^Mp, = IQ-^-^m^ and zx = Mp/cx and = zx/2. Results 
for (cx,c<^) = (30,2) and (40,5) are shown together with the analytic estimate based on 



Next, we consider the other case : Hini ~ as in the chaotic inflation model in 
supergravity ^22j. The superpotential is given by 

W = m^Xcf), (36) 

and the inflaton starts to oscillate from about the Planck scale. The X is stabilized at 
the origin during and after inflation. We therefore take 

- %■ 

In this case, the modulus oscillates in a potential which changes with a frequency m'^^. 
Since the modulus mass and the inflaton mass become comparable for a certain time, 
it is likely that resonant particle production occurs, making it difficult to estimate the 
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modulus abundance analytically. Let us therefore give a very rough estimate on the 
abundance of the moduli in the form of coherent oscillations. We assume cx = 0(10) and 
= 0{1). Since both Hi and H2 oscillates at a frequency of with an opposite phase, 
the modulus mass becomes of order the Hubble parameter when Hi vanishes. During the 
time of H~^ ~ m^, the modulus moves from zx to ~ z^. As the modulus mass becomes 
greater than the Hubble parameter, it comes back to zx again and starts oscillating with 
an initial amplitude given by ~ \Z(i, — zx\- Therefore, the Polonyi abundance in the form 
of coherent oscillations is estimated as 

^-ItJ^^C^] (39) 



for Cx = 0(10) and = 0{1). Numerical results are shown in Fig. El Here the modulus 
abundance is evaluated as a function of cx in the multi chaotic inflation model. We 
have taken m,^ = O.lMp, = 10~^Mp, = lO'^'^m^, Ccf, = 2 and zx = Mp/cx and 
= zxj^. One can see that there is a resonant feature and the analytic estimate roughly 
agrees with the numerical result. Note however that we have focused on the homogeneous 
mode of the modulus, and the result might be significantly affected by taking account 
of the resonant particle production, or preheating. Assuming that the resonant particle 
production, if any, is just added to the above estimate, the modulus abundance is huge 
and another solution to the cosmological Polonyi/moduli problem is needed,^ 

After all, in the multi-field inflation models, the adiabatic suppression is inefficient 
once the modulus production at the end of inflation is taken into account. We need to 
tune the couplings of the multi-inflaton fields to the Polonyi/moduli so that zx — in 
order to solve the cosmological Polonyi/moduli problem, although the required amount 
of tuning is relaxed from 10"^'^ to 10~^ compared with the case of c = 0{1). 

3.4 Required fine-tuning to solve the Polonyi/moduli problem 

In order to see the required amount of tuning more quantitatively, we define the tuning 
factor e as \zx — z^\ = eMp and examine how small e must be in order to solve the 
cosmological moduli problem. To be concrete, we rely on the estimate ( 132|) . This estimate 



If C0 is much smaller than unity for some reason, the modulus abundance is likely suppressed. 
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Figure 6: The modulus abundance as a function of cx in the muhi chaotic inflation model. 



We have taken = O.lMp, 
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-2.5 



= 2 and zx = Mp/cx 



and = zx/2. Analytic estimate based on (139|) is also shown. 



applies almost identically to the case of single- field inflation models if we set c,^ = 1, as 
shown in Eq. ( 122|) . Fig. [7] shows the contours of the required amount of the fine-tuning 
e for the adiabatic suppression to solve the moduli problem, on the plane of rriz and Tr. 
Here we have taken ifjnf = 10^° GeV and cx = 30, = 1. The required amount of tuning 
becomes less (more) stringent by one order of magnitude for Hi^f = 10^^ (10^) GeV. 
That is to say, e scales as (ifinf/10^°GeV)^/^. In the shaded region at the upper left 
corner, the adiabatic suppression does not work. We have used constraints from BBN, 
diffuse X(7)-ray and overclosure as done in Ref. [25] with updated BBN constraints. For 
rriz > 100 TeV, we have assumed anomaly mediation relation between the gravitino and 
the Wino-LSP mass and derived constraints from LSP overproduction by the modulus 
decay. We assumed that the modulus has couplings with the SM particles as Eq. (153|1 
with C3 = 1. We have also assumed that the modulus mass is comparable to the gravitino 
mass so that the gravitino production from the modulus decay is inefficient. Notice that 
upper bounds on the reheating temperature coming from Polonyi/moduli and gravitino 
thermal production, which will be discussed in the following sections, are not taken into 
account in this figure. From this figure, it is seen that for wide range of modulus mass and 
reheating temperature, the required amount of tuning lies in the range of 10~^ ^ e ^ 1- 
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[GeV] 

Figure 7: Contours of the required amount of the fine-tuning e for the adiabatic sup- 
pression to solve the moduh problem, on the plane of niz and Tr. Here we have taken 
-f^inf = 10^° GeV and cx — SO, c^f, — 1. The required amount of tuning becomes less 
(more) stringent by one order of magnitude for Hi^f — 10^^ (10^) GeV. That is to say, e 
scales as {Hi^f/lO^^GeVy/^. In the shaded region at the upper left corner, the adiabatic 
suppression does not work. Notice that upper bounds on the reheating temperature com- 
ing from Polonyi/moduli and gravitino thermal production, which will be discussed in the 
following sections, are not taken into account in this figure. 
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3.5 Other form of the modulus potential 

Before closing this section, we comment on whether the adiabatic suppression works in 
other types of the modulus potential. First we consider the case of a negative Hubble 
mass squared. The scalar potential is 

V{z) = -{ml-c^H^)z^ + ^^z'', (40) 

where n{> 4) is an even integer and A is a positive constant. The flat directions in SSM 
have the potential of this type [5]. The true minimum is obviously z = Q. The temporal 
minimum is given by 

.„„(i/)=|f?"<^'"'-'"5'/^l'"""' t^"' . (41) 

1^ for ciJ < m^. 

Hence we obtain 

f{H)= ^ ^ for ci7>m„ (42) 

in the matter-dominated phase. Therefore, f{H) diverges aX H = m^jc and the adiabatic- 
ity is always violated independently of the size of c and A. Thus the adiabatic suppression 
does not occur for the scalar field of the potential of (HOl) . The z field begins to oscillate 
at ~ ^z/c around the minimum z = with amplitude of ~ (M^-^m^/A)!/^"-^). We 
have numerically confirmed that the oscillation amplitude is not suppressed regardless of 
the value of c. 

The situation slightly changes if z has the negative mass term around the origin : 

V{z) = -{-ml-c^H^)z^ + ^^z''. (43) 

The temporal minimum is given by 

Zm,.{H) = [M--\c'H' + ml)/X] ^/^""^^ . (44) 

It smoothly connects to the true minimum Zj^^^ = (M^-^m2/A)'/^""'^ In this case we 
obtain 

f(H) = \ ^ (45) 

in the matter-dominated phase, hence the time scale of the potential change is roughly 
given by H. This is always smaller than the mass scale around the temporal minimum 



25 



for c ^ 1 and hence the adiabaticity condition ( JTTl) is met in the oscillation regime. As is 
already explained, however, z cannot track the temporal minimum at the very beginning 
of its motion. Thus the coherent oscillation is necessarily induced as explained before. 
The induced modulus amplitude is estimated as 



-1 



(46) 



and the resulting modulus abundance is 



Pz_ f ZmmjHi) 

s ^ 8 Mp 



(47) 



This is the minimum modulus abundance in this model. By tuning the initial velocity 
i, we can suppress the abundance further, but such a tuning is not likely to occur in a 
realistic setup starting from inflation. Numerical calculations and analytic estimate are 
shown in Fig. [HI We have taken n = 6 and A = Mp/ml = 10^ and Hi = lO^m^. Analytic 
estimate does not depend on c for n = 6. Our estimate fits very well with the numerical 
results. 
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Figure 8: Time evolution of the modulus abundance for c = 1, 10, 100 for the model fl43|) . 
We have taken n = 6 and A = Mp/ml = 10^, = 10~^-^m^ and Hi = lO^m^. 



The last example is a very flat potential such as 

V = Vo log (1 + ^) + Ic'H^z - z^f. (48) 
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Time evolution of the modulus abundance for c = 1, 10 for the model (HHl) is shown in 
Fig. H We have taken = Mp/c, M = zJlO and Vq = rri^zl with m = lO^^Mp. 
Adiabatic suppression does not work in this case, simply because the mass around the 
temporal minimum changes the sign at some z < z^ and the adiabaticity is necessarily 
violated. If, instead, we have chosen M ^ z^, the effective mass is always positive for 
z < z^ and the model is effectively described by ([1]), hence the adiabatic suppression 
occurs for c ^ 1. 




t 

Figure 9: Time evolution of the modulus abundance for c = 1, 10 for the model fHHl) . We 
have taken z^ = Mp/c, M = z^/10 and Vq = m?zl with m = 10"^Mp. 



3.6 Summary 

To summarize, it depends on the behavior of the time-dependent potential minimum 
whether the adiabatic suppression takes place. In particular, if the curvature of the 
potential about the minimum vanishes temporarily, or if the position of the minimum 
changes rapidly compared to the curvature, the adiabatic solution does not work. There- 
fore it should be kept in mind that the adiabatic solution does not apply to any moduli 
with an arbitrary potential. Nevertheless, it is remarkable that the adiabatic suppression 
mechanism works for the simplest example like ([1]), which includes the Polonyi field. 
It should be noted, however, that a considerable amount of modulus oscillation is 
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generically induced at the end of inflation, as already shown. Some amount of tuning is 
needed in order to solve the cosmological moduli problem for relatively high reheating 
temperature, although the required amount of tuning is significantly reduced compared 
with the ordinary case without the adiabatic suppression. 

In the following sections we always assume that the modulus/Polonyi coherent oscil- 
lation is somehow suppressed by the adiabatic suppression and some amount of tuning if 
needed and discuss how the reheating temperature is constrained even in such a case. 

4 The Polonyi problem in gravity mediation 

Let us consider the Polonyi problem in gravity mediation, where SUSY is broken by the 
F-term of the Polonyi field Z . We consider the following Kahler potential, 

+ (49) 

where represents the potential minimum during inflation. Here we have omitted in- 
teractions like K D nMpZ + K,'\(f)\'^ZZ/Mp + ■ ■ ■. The linear term is necessary to set 
the origin to be the low-energy minimum, which also shifts the minimum during inflation 
from by some factor, but not an order of magnitude. So we drop the linear term since 
it does not change the modulus dynamics significantly in our context. The latter coupling 
induces the inflaton decay into a pair of the gravitinos [27], which may impose severe 
constraints on inflation models if ^ 1. 

The coupling with the inflaton in the Kahler potential (H^ induces a Hubble mass 
correction to the effective potential of the Polonyi field as 

V{z) D c^H^\z-z,\'^ (50) 

where z denotes a scalar component of Z, with 

c' = 3(c? + l), (51) 

which is equivalent to Eq. ([T]). The Hubble-induced mass term is present during inflation 
and the inflaton oscillation era, but it disappears after the reheating. The Polonyi mass 
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around the origin in the low energy is 







ml = (3c2 + 2y3)m2/2, (52) 

In order to suppress the coherent oscillations of the Polonyi field, we assume c = 
C(10) and Tr ^ 0.05a/?t7^Mp. These two conditions must be satisfied for the adiabatic 
suppression to work, independently of the details of the infiaton dynamics. As shown in 
the previous section, the Polonyi coherent oscillation is generically induced at the end of 
infiation. Hereafter we simply assume that this contribution is suppressed at the price of 
fine-tuning of at most 10""^ level (see Fig. [7]). 

The Polonyi couplings to the SSM fields are parametrized as 

7 + cai^) WaW'' + h.c, (53) 



4 V "M, 



i^e{-^}^J^_Z\f? + i....]\. (54) 



Mf, \Mp 

where Wa is a field strength of the gauge supermultiplet, / collectively denotes the chiral 
matter superfield, and C3, C4 and C5 are numerical coefficients. In the following we take 
all the coupling constants Ci, ■ ■ ■ , C5 to be real, for simplicity. These terms generate the 
soft SUSY breaking masses as 

m-g = m3/2, (55) 

m'j = {cl + 4 + l)ml/2, (56) 

where rrig and denote the gaugino and sfermion masses, respectively. The above 
couplings enable the Polonyi field to decay into the SSM particles, and more important, 
the Polonyi is necessarily produced from particle scattering in thermal plasma through 
the couplings with C3 and C5. (Note that the origin of Z is set to be the low-energy 
potential minimum.) Therefore, even if the coherent oscillations of the Polonyi field can 
be negligibly small by the adiabatic suppression mechanism, its abundance may be still 
non-negligible if the reheating temperature is high. 



Precisely speaking, this is the mass of the real component of Z in the presence of the linear term 
with K = \/3 — 1. The imaginary part has a slightly smaller but comparable mass. 
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parameter 


definition 


effects 


Cl 


Eq. (HHD 


Hubble mass for tlie Polonyi 


C2 


Eq. (HHD 


Polonyi mass 


C3 


Eq. dSSD 


Gaugino mass, Polonyi abundance 


C4 


Eq. dMD 


Sfermion mass 


C5 


Eq. dSl 


Sfermion mass, Polonyi abundance 



Table 1: Coefficients of non-renormalizable operators and their effects on masses and the 
Polonyi abundance from thermal scattering. 

In the rest of this section, we consider cases with several different values of Cj. First we 
consider the case that only ci is enhanced while all the other couplings are of order unity. 
This is the minimal set-up to solve the Polonyi problem using the adiabatic solution. Next 
we extend the minimal set-up to allow enhancement of other couplings. Such extension 
may be indeed reasonable; if the enhanced coupling of Ci is due to an exchange of fields 
of mass below the Planck scale or due to some strong dynamics at the Planck scale, we 
naively expect a coupling like C2 is similarly enhanced. Furthermore, if the SSM particles 
are involved with the strong coupling at the Planck scale, or if the fundamental cut-off 
scale of theory is one order of magnitude smaller than the Planck scale, we expect that 
all the modulus couplings to the SSM particles, C3, C4 and C5 as well are enhanced. As we 
will see below, the constraints on the reheating temperature depends on the size of these 
couplings constants. 

4.1 The minimal set-up 

First we consider a case that the Polonyi-inflaton coupling (ci) is enhanced to solve the 
Polonyi problem, while all the other couplings among the Polonyi and SSM particles are 
assumed to be Planck-suppressed with coefficients of order unity. 

The interaction f lS^ induces the Polonyi decay into the SM gauge bosons, which is 
subject to stringent constraints from BBN [T71 [281 129] . The partial decay rate is estimated 
as 

where we have considered the decay into all the SM gauge bosons. The partial decay rates 
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^Z^gg - ^3 32^^2 ("^3/2 + , (58) 



of the Polonyi into SSM gauginos is given by, if kinematically allowed, 

'p 

where we have approximated all the gaugino masses are same. The partial decay rates of 
the Polonyi into sfermions is also close to the above rate. Note that, if the Polonyi mass 
is much heavier than the gravitino mass (i.e. C2 ^ 1), the decay into the gauge bosons is 
the dominant decay mode [5U] . 

The Polonyi is necessarily produced by thermal scattering like the gravitino, even if 
the coherent oscillations are suppressed by the adiabatic suppression mechanism. The 
gravitino abundance is [311 132] 

---=f— "(-4;)(io4v)' (-) 

where 113/2 is the number density of the gravitino, and s the entropy density. The Polonyi 
abundance is expected to be the same order of that of the transverse component of the 
gravitinos. So we use the following estimate 

.r.-l2('2 , „2n ( 



y,.^.10-(e^,e|)^^j. (60) 

where nz is the number density of the moduli. Note that there is an uncertainty of order 
unity in the above estimate. 

Using the Polonyi and gravitino abundances, we can derive cosmological constraints 
on the reheating temperature. We show in Fig. [10] upper bounds on the reheating temper- 
ature as a function of the gravitino mass. Here we set = 7/13/2, C3 = 1(5) and C5 = 1 in 
the top (bottom) panel. The precise value of C4 is not relevant for the following arguments 
as long as \ci\ < 0{1). Here and in what follows the GUT relation among gaugino masses 
is assumed, unless otherwise stated. In the top (bottom) panel, the LSP is the bino- or 
higgsino-like neutralino (gravitino). The constraints on the thermal relic abundance of the 
SSM LSP are not taken into account in the both panels, which will be discussed shortly. 
The meaning of each hue is as follows. "BBN (TH Polonyi)" and "CMB (TH Polonyi)" 
refers to the BBN [T71 |28l [29] and CMB [33] bounds on the thermally produced Polonyi, 
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"LSP from TH Polonyi" to the LSP overproduction bound from the Polonyi decay|j "TH 
Gravitino" to the bound from the gravitino thermal production, taking account of the 
gravitino decay effects on BBN and LSP overproduction, and "Adiabaticity" to the bound 
for successful adiabatic suppression on the coherent oscillation of the Polonyi. We can 
see from Fig. [10] that the Polonyi problem can be solved for a sufficiently low reheating 
temperature, Tr < 10^ GeV. The adiabatic solution is therefore an attractive solution to 
the Polonyi problem. On the other hand, the reheating temperature above 10^ GeV is not 
allowed, and therefore the thermal leptogenesis does not work [10], if all the constraints 
are taken at face value. 

Now let us discuss if some of the constraints could be relaxed. In the case of the 
bino/higgsino LSP (top panel), the BBN constraint on the gravitino decay is so stringent 
for m3/2 ^ lOTeV. So let us focus on the heavy gravitino of mass > 10 TeV, where the 
most stringent constraint comes from the LSP overabundance (dotted pink). In order to 
weaken the bound, we need to suppress both the thermal and non-thermal relic of the 
LSP. (Note that the thermal relic abundance of the LSP is not taken into account in 
the figure.) In order to suppress the thermal relic abundance, the higgsino component 
should be sizable. The non-thermal component can be suppressed if the LSP mass is 
much lighter than the gravitino. Thus, one possibility to relax the constraint in the heavy 
gravitino region is to assume that the LSP is the light mixed bino-higgisino neutralino, 
namely, mbino ~ ^ "^3/2, which requires C3 ^ 1. Then the SUSY mass spectrum would 
be similar to the focus point region [31], which is usually considered as a region where 
thermal leptogenesis works when only the gravitino constraint is considered. The bound 
can be similarly relaxed for the (purely) higgsino-like LSP, much lighter than the SSM 
gauginos and sfermions. In both cases the bound can be relaxed by an order of several 
tens to hundred, and the reheating temperature can be as high as O{10^) GeV. Another 
possibility to relax the bound is to introduce R-parity violation; then the constraint 
disappears, but another DM candidate such as an axion would be needed. 

In the case of the gravitino LSP (bottom panel), the BBN constraint on the Polonyi 

* Note that in this figure we do not take into account the thermal relic abundance of the LSP. In order 
to avoid the overabundance, the LSP mass cannot be arbitrarily large, although the precise thermal 
abundance depends on other SUSY particle masses. 



32 



Tr [GeV] 

-,12 



Polonyi : = 1,C3 = 1 , Bino LSP 




BBN (TH Polonyi) " 

CMB (TH Polonyi) 

LSP from 
TH Polonyi 

Adiabaticity 
TH Gravitino ■ 



m3/2 [GeV] 



Tr [GeV] 

10^2, 



Polonyi : = \,c^ = 5 , Gravitino LSP 



10 
10^ 
10^ 
10^ 
10^ 
10° 



10 



BBN (TH Polonyi) ■ 

CMB (TH Polonyi) 

LSP from 
TH Polonyi 

Adiabaticity 
TH Gravitino 



10' 



10"^ 



lO"" 



m3/2 [GeV] 



Figure 10: Upper bounds on the reheating temperature for the Polonyi model with 
T^z = "^3/2 (c2 = 1) for C3 = 1 (top) and C3 = 5 (bottom). In the top (bottom) panel, 
LSP is the bino/higgsino (gravitino). "BBN (TH Polonyi)" and "CMB (TH Polonyi)" 
refers to the BBN and CMB bounds on the thermally produced Polonyi, "LSP from TH 
Polonyi" to the LSP overproduction bound from the Polonyi decay, "TH Gravitino" to 
the bound from the gravitino thermal production, taking account of the gravitino decay 
effects on BBN and LSP overproduction, and "Adiabaticity" to the bound for successful 
adiabatic suppression on the coherent oscillation of the Polonyi. Note that in this figure 
we do not take into account the thermal relic abundance of the bino/higgsino. 
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decay is so stringent that it is difficult to make the thermal leptogenesis viable. This 
should be contrasted to the usual case in which the reheating temperature can be as high 
as C(109) GeV at mg/s = O{10^) GeV in the absence of the Polonyi field. 

In summary, the adiabatic solution, with fine-tuning on the Polonyi potential at 0.1% 
level, works in the minimal set-up if Tr is sufficiently low, Tr < 10^ GeV. The constraint 
on the reheating temperature can be relaxed in several ways. In the case of the neu- 
tralino LSP, the thermal letpogenesis becomes possible for the gravitino mass as heavy 
as (9(10) TeV, if the lightest neutralino is the mixed bino-higgsino or purely higgsino-like 
LSP, or if the R-parity is violated. In the case of the mixed bino-higgsino LSP, the SUSY 
mass spectrum is similar to the focus point region. On the other hand, the thermal lepto- 
genesis [10] does not work in the case of the gravitino LSP because of the stringent BBN 
bound on the Polonyi decay. 

4.2 Extended adiabatic solution 

Now let us extend the minimal set-up to allow other parameters to be enhanced. In 
particular we focus on the case of Ci ^ 1 and C2 ^ 1. Such an enhancement of C2 is 
indeed expected in a certain theoretical framework. In this case the Polonyi is much 
heavier than the gravitino and the cosmological problems associated with the thermally 
produced Polonyi can be relaxed. 

In addition to ci and C2, the Polonyi may also couple to SSM fields with enhanced 
interactions. Actually, as shown in Ref . [12] , it is conceivable that all the Polonyi couplings 
are universally enhanced, ci, C2, C3, C4, C5 = C(10). The gaugino mass is determined by 
C3, while the sfermion mass depends on C4 and C5. If C3, C4, and/or C5 are larger than 
0{1), the gauginos and sfermions are heavier than the gravitino. On the other hand, the 
Polonyi abundance (160|) depends on C3 and C5. In the following we vary those parameters 
and analyze how they affect the cosmological bounds. 

The Polonyi is much heavier than the gravitino, if C2 ^ 1. Then the Polonyi decays 
into a pair of the gravitinos, as well as into the gauge bosons^ The decay rate into the 
^ Decay into gauginos are less efficient for 7713/2 <C rriz [3D1. 
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gravitino pair is given by 

= ^^;;t^- (61) 

In fact, Z — )■ 2'?/'3/2 is often the dominant decay mode. If this is the case, the branching 
fraction of decay into visible particles is suppressed, which relaxes the BBN constraint on 
the Polonyi decay if the gravitino is the LSP and stable. 

The Polonyi and gravitino abundances, given by Eqs. fIBU]) and f lS^ . respectively, are 
enhanced for C3 ^ 1, namely, rrig 3> Although the partial decay rate of the Polonyi 

field into gauge bosons is enhanced for C3 ^ 1, the branching fraction remains small if 
the gravitino production is the main decay mode. 

The upper bounds on the reheating temperature are shown for C2 = 30 and C3 = 
C5 = 1(5) in the top(bottom) panel in Fig. [11] The meaning of each line is same as in 
Fig. [ini In the top panel, the LSP is either bino- or higgsino-like (or mixed bino-higgsino) 
neutralino, while the gravitino is the LSP in the bottom panel. Compared to Fig. [TOl 
the gravitino constraint is same as before in both cases, while the BBN constraint on the 
Polonyi decay is significantly relaxed. This is because the Polonyi mass is much heavier 
than the gravitino and it decays before BBN begins, or at an early stage of BBN where 
the constraint from the helium overproduction is relatively weak. 

In the top panel, the overall constraints are similar to Fig. [TDl and the constraint 
from the LSP overabundance in the heavy gravitino region can be relaxed in a similar 
way to the previous case. The reheating temperature as high as 10^ GeV is allowed 
only for the heavy gravitino case of ^3/2 ~ lOTeV. Even if we demand a moderate 
reheating temperature of Tr > 10^ GeV for non-thermal leptogenesis [HI |35], we need 
a relatively heavy gravitino, m3/2 ^ 5 TeV. On the other hand, in the bottom panel, 
since the BBN constraint on the Polonyi field is significantly relaxed, there appears an 
interesting parameter region at ^3/2 ~ 100 GeV with Tr ~ 10^ GeV where the thermal 
leptogenesis may work successfully. Here it should be noted that C3 cannot be much larger 
since it enhances the gaugino mass and hence the gravitino thermal production. If we 
demand Tr > 10^ GeV for non-thermal leptogenesis, a broader region for the gravitno 
mass (9(10)GeV-(9(10)TeV is allowed. In this region gravitino is the LSP, while sfermions 
are heavier than the gauginos. This resembles the mass spectrum in the so-called focus- 
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Figure 11: Same as Fig. [10] but for C2 = 30 and C3 = 1 (top) and C2 = 30 and C3 
(bottom) with C3 = C5. In the top (bottom) panel, LSP is the bino (gravitino). 



Note that the effect of the next-to-hghtest SUSY particle (NLSP) is not taken into 
account in the figure. In fact, the long-lived NLSP significantly affects BBN in the case 
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of gravitino LSP and hence we need a small i?-parity violation in order for the NLSP to 
decay well before BBN begins. The NLSP lifetime depends on the pattern of i?-parity 
violation and its magnitude. In particular, if the i?-parity is violated by the tri-linear 
interaction such as u'^d^d'^ or LLe^, the NLSP lifetime depends crucially on the sfermion 
masses. The displaced vertex may give useful information on discriminating the model. 

5 The moduli problem in gauge-mediation 

The gauge-mediated SUSY breaking (GMSB) models have attractive features that 
they do not suffer from the SUSY flavor problem. In the GMSB model, the Polonyi- 
like field is not always required in the SUSY breaking sector. Thus there may not be a 
Polonyi problem. But even in this case there may exist modulus fields especially if the 
four dimensional supergravity arises from the compactification of the extra dimensions in 
string theory. These moduli are expected to have masses of the gravitino and Planck- 
suppressed interactions, and hence they cause serious cosmological problems. 

In this case the relevant parameters are the infiaton-modulus coupling ci, which must 
be enhanced in order to solve the moduli problem. As shown in Fig. [TJ the moduli 
oscillation can be sufficiently suppressed by the adiabatic solution without fine-tuning for 
niz < 1 MeV. However, we still need to care about thermally produced moduli. In the 
GMSB model, the moduli, as well as the gravitino, are light and the lifetime are longer 
than the case of gravity-mediation. If the modulus survives after the recombination epoch, 
the decay produced photons contribute to the diffuse X(7)-ray background [37] which may 
easily exceed the observational limit [381 ES] • 

Fig. [T2l shows constraints on the reheating temperature in the GMSB model with mod- 
uli. We take = (top) and = 10m3/2 (bottom) and rng=l TeV. "BBN (TH 
moduh)", "CMB (TH moduh)" and "Diffuse gamma (TH moduh)" refer to the BBN, 
CMB and diffuse gamma-ray flux bounds on the thermally produced moduli, "Adiabatic- 
ity" to the bound in order for suppress the Polonyi abundance and "TH Gravitino" to 
the usual bound from the gravitino thermal production. 

First it is remarkable that the ultra-light gravitino mass region, 7713/2 < 16 eV, where 
it was believed that there is no upper bound on the reheating temperature, is clearly 
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inconsistent with high-reheating temperature scenario once we demand that the adiabatic 
mechanism suppresses the modulus abundance!!^ 

For the case of intermediate gravitino mass of (9(10) GeV, the situation is similar to 
that studied in the previous section. It may be consistent with non-thermal leptogenesis 
scenario if the modulus mass is enhanced compared to the gravitino mass under the broken 
i?-parity. 

6 The moduli problem in anomaly-mediation 

Finally we mention the moduli problem in the anomaly-mediated SUSY breaking (AMSB) 
model [3]. As noted in Introduction, the moduli problem is milder than that in the case 
of gravity- or gauge-mediation since the moduli are heavy enough to decay before BBN. 
But it was recognized in Ref. [1] that the gravitinos produced by the modulus decay cause 
cosmological problems. Even if the gravitinos are also heavy enough, the abundance of 
LSPs produced by the modulus/gravitino decay is too much. A possible solution is to 
introduce a small i?-parity violation in order for the LSP to decay. In any case, the 
modulus decay produces huge amount of entropy and the preexisting baryon asymmetry 
is significantly diluted!"] 

Here we consider the adiabatic solution to the moduli problem in the AMSB model, 
which is free from the problems of late-time entropy production and hence the dilution of 
the baryon asymmetry. As shown in Fig. [71 we need a tuning on the modulus potential 
at the level of O{10^^) or (9(10~^) in order to solve the moduli problem in the original 
sense for a relatively high-reheating temperature consistent with thermal or nonthermal 
leptogenesis. Again, however, we still need to take account of thermal production of the 
moduli. Fig. [T^ shows upper bounds on the reheating temperature in AMSB models with 
moduli. We take rriz = rn3/2 (top) and = 10m3/2 (bottom). The meaning of each 
line is same as those in Fig. [121 except for the "LSP from TH moduli" which denotes the 

The smallness of the modulus abundance in this mass range may be explained by the anthropic 
arguments without invoking the adiabatic suppression, since in this region the constraint comes from the 
modulus overabundance as the DM. 

The amount of dilution factor is roughly given by ~ Tz/T^ where Tz is the modulus decay tempera- 
ture. The most baryogenesis models do not work under this condition, except for some parameter ranges 
in the Affleck-Dine mechanism [6]. 
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Figure 12: Upper bounds on the reheating temperature in GMSB models with moduh. 
We take ruz = (top) and = 10m3/2 (bottom) and Tng=l TeV. "BBN (TH mod- 
uh)", "CMB (TH moduh)" and "Diffuse gamma (TH moduh)" refer to the BBN, CMB 
and diffuse gamma-ray flux bounds on the thermally produced moduli, "Adiabaticity" to 
the bound in order for suppress the Polonyi abundance, and "TH Gravitino" to the usual 
bound from the gravitino thermal production. 
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LSP overproduction bound from the decay of thermally produced moduli. Note that in 
this case the "TH gravitino" bound includes the LSP overproduction from the gravitino 
decay. It is found that the gravitino mass of around lOOTeV is compatible with the 
thermal leptogenesis scenario (Tr > lO^GeV). We do not need an i?-parity violation in 
this case. 

Notice that the Wino is the LSP in AMSB models and its mass is given by = 
((71/1671^)7713/2 ~ 2.6 X 10"'^m3/2- In the parameter region considered above, the Wino 
produced non-thermally by the gravitino/modulus decay can be the dominant component 
of DM. But the Wino has rather large annihilation cross section and the Wino DM mass 
is limited from the observation of light element abundances [lOl HI], cosmic microwave 
background anisotropy [42l |43] and gamma-rays by the Fermi satellite |44] . These con- 
straints demand m^^ > 200 GeV which translates into the bound on the gravitino mass 
as 777.3/2 ^ 80 TeV if the DM mainly consists of the Wino. 

7 Conclusions and discussion 

One of the major obstacles to construct a consistent cosmological scenarios in the most 
SUSY breaking models is the cosmological Polonyi/moduli problem. We have examined 
carefully the adiabatic solution to the cosmological Polonyi/moduli problem, and found 
that the Polonyi/moduli oscillation is necessarily induced at the end of inflation, and its 
abundance depends on the inflation models. As a result, the Polonyi/moduli problem 
cannot be solved without further tuning on the potential parameters for a broad range 
of the modulus mass. The possible ranges of the moduli mass where no severe tuning 
is needed for solving the moduli problem are either rriz < 1 MeV or rriz > 10 TeV (see 
Fig. [7j). In the former case, the reheating temperature cannot be larger than ~ ITeV 
from the gravitino constraint (see Fig. [T2|) . In the latter case, the reheating temperature 
can be as high as ~ 10® GeV with a mild tuning of (9(0.1 — 0.01) and may be consistent 
with non-thermal leptogenesis scenario [211 ES]- If we allow a tuning at the 1% level, 
the thermal leptogenesis requires rriz ^ 100 TeV. On the other hand, for a very low 
reheating temperature of (9(1) MeV, the adiabatic suppression mechanism can solve the 
Polonyi/moduli problem without tuning for every range of Polonyi /modulus mass. 
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Figure 13: Upper bounds on the reheating temperature in AMSB models with moduh. We 
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We have also pointed out that Polonyi/moduh production from scatterings in thermal 
bath can be relevant, even if the coherent oscillations are suppressed by the adiabatic 
mechanism. If the Polonyi/moduli have Planck-suppressed interactions with the SSM 
particles, their abundance is comparable to that of the transverse component of the grav- 
itino. In fact, the cosmo logical constraint on the Polonyi field produced from thermal 
scattering has turned out to be stringent for a broad range of the gravitino mass (see 

Fig.nni). 

Since the adiabatic solution necessitates the enhanced Polonyi/moduli- inflaton cou- 
pling, other couplings may also be enhanced, depending on the origin of the enhanced 
coupling. We have studied several cases where the Polonyi self-coupling (02) and/or the 
Polonyi-matter couplings (03, C4, C5) are enhanced. In the gravity-mediation model, the 
Polonyi field mainly decays into a pair of gravitinos, if its self-coupling is enhanced. This 
relaxes the cosmological constraint on the Polonyi abundance when the gravitino is the 
LSP. Once we allow a fine-tuning on the Polonyi potential at the level of O{10~^) or 
O{10~^) , the adiabatic solution works in consistent with thermal or nonthermal leptoge- 
nesis scenario, for some gravitino mass ranges (see Fig. [7]). In particular there appears a 
new interesting parameter space where the gauginos are relatively light (~ (9(100) GeV) 
and the gravitino is the LSP. A small amount of R-parity violation is needed to avoid 
the BBN constraint on the NLSP decay. Sfermion masses can be comparable to or much 
heavier than the gaugino mass depending on the parameters C4 and C5. We have also com- 
mented that non-thermal gravitino production from the inflaton decay may be enhanced 
in this setup. 

In the gauge-mediated SUSY breaking model, there may also exist light moduli, and 
we have studied the adiabatic solution to the moduli problem. It is found that the ultra- 
light gravitino scenario {1713/2 < 16eV), where the gravitino problem does not exist, is not 
consistent with the adiabatic solution unless Tr < 1 TeV, although the modulus amplitude 
may be anthropically tuned without the adiabatic suppression mechanism. 

In the anomaly-mediation model, once the modulus oscillation is suppressed by the 
adiabatic suppression mechanism and some tuning, constraints from thermally produced 
moduli are rather weak. The reheating temperature can be high enough to be consistent 
with (non-)thermal leptogenesis. 
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We conclude that there are many non-trivial constraints for the adiabatic suppression 
on the moduli to work successfully. Even in the presence of adiabatic suppression, we 
need a tuning on the modulus potential in order to avoid the moduli problem especially 
for a relatively high reheating temperature (see Fig. [7j), although the required amount of 
tuning is significantly relaxed. For a low reheating temperature, on the other hand, the 
adiabatic suppression mechanism can solve the Polonyi/moduli problem without any fine- 
tuning. Thermal production of the moduli also set a severe constraint on the reheating 
temperature for relatively light moduli of < 1 TeV. These aspects of the moduli 
should be taken into account when the solution to the moduli problem in the adiabatic 
suppression is discussed. 
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